How edge states are destroyed in disordered mesoscopic samples? 
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We report theoretical investigations on how edge states are destroyed in disordered mesoscopic 
samples by calculating a "phase diagram" in terms of energy versus disorder strength {E,W), and 
magnetic field versus disorder strength {B,W), in the integer quantum Hall regime. It is found 
that as the disorder strength W increases, edge states are destroyed one by one if transmission 
eigen-channels are used to characterize the edge states. Near the insulating regime, transmission 
eigen-channels are closed one by one in the same order as edges states are destroyed. To identify 
those edge states which have survived disorder, we introduce a generalized current density that can 
be calculated and visualized. 
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When a two-dimensional mesoscopic sample is sub- 
jected to external magnetic field, peculiar electronic 
states — edge states, may be established at the bound- 
aries of the sample pj, i2] . Classically, Lorentz force pushes 
electrons toward the sample boundary and electron tra- 
jectories become skipping orbits. Edge states can be con- 
sidered as the quantum version of skipping orbits 2]. Im- 
portantly, edge states in mesoscopic samples provide nec- 
essary density of states (DOS) between the Landau levels, 
integer quantum Hall effect (IQHE) can therefore occur 
in the clean sample limit 2]. In contrast, for infinitely 
large samples, a degree of disorder in the sample appears 
necessary which provides DOS in between Landau levels 
to stablize the Fermi energy for IQHE^. Nevertheless, 
increasing disorder will eventually destroy IQHE and how 
does this happen has been an important issue attracting 
numerous studies. 

Here we address the disorder issue for mesoscopic sam- 
ples, namely how edge states are destroyed by disorder in 
the IQHE regime and, at a fixed filling factor v > 1, are 
edge states destroyed all at once or one by one. These 
important questions provide insight to the IQHE phase 
diagram for mesoscopic samples, and may shed hght to 
similar problems in samples of infinite size. We address 
these questions by extensive calculations on a mesoscopic 
graphene system and a square lattice model (see inset of 
FigUJi') to map out a "phase diagram" of edge states in 
the presence of disorder. Here the "phases" in the "phase 
diagram" denote quantum states and no phase transi- 
tions are implied between these states. We use transmis- 
sion eigen-channels to characterize edge states (see be- 
low), and we found that they are destroyed one by one. 
At large enough disorder, the system becomes an insu- 
lator and transmission eigen-channels are closed one by 
one in the same order as the edge states are destroyed. 

We begin by discussing our definition of edge states 
as well as the way to visualize them. In transport the- 
ory, for each incoming channel of a semi-infinite lead a 
whose wave function is \Wa,m > where m = 1, 2, • • • iV 
denotes one of the N channels, one solves a scattering 



problem. iW^a,™ > is an eigen-channel of lead-a but not 
the entire device. To find the eigen-channels of the en- 
tire device, we diagonalizejj] the transmission matrix T 
by a unitary transformation U, for a two probe device 
having scattering matrix Srl, T^ — Sj^j^Srl- Mathe- 
matically, this means acting U on the incoming chan- 
nels \Wa^m > (which is a column vector with N compo- 
nents) to obtain a new set of orthogonal incoming modes 

\Wam >== Y.n \^a,n > Unm- OnCe doue, \Wa,m > IS 

an eigen-state of the entire device (leads plus scattering 
region). In other words, if an incoming electron comes at 
state \Wa,m >, it will traverse the entire device without 
mixing with any other eigen-state \Wa,m' >■ This way, 
the resulting transmission matrix T = WTU is diago- 
nal. In the presence of a strong magnetic field B, it is 
therefore natural to identify \Wa,n{B) > as edge states 
because they are the eigen-states or eigen-channels of the 
entire device sample. How to visualize edge states in the 
IQHE regime? This may be achieved by plotting the 
current density in real space. The subtle issues of cur- 
rent density in IQHE has been discussed in Refja. In our 
work where there are disorder in the sample, the eigen- 
channels provide a convenient way to define a general- 
ized current density for each channel — since the eigen- 
channels do not mix. Clearly, the total transport cur- 
rent is obtained by integrating current density along any 
cross-section perpendicular to the transport direction. 

The total transmission coefficient is obtained from the 
N X N Hermitian transmission matrix T = {ty}. Ap- 
plying an unitary transformation C/[J|, we obtain T = 
U^TU — {tiSij}, which is diagonal with elements ti. The 
transmission coefficient ti of eigen-channel i is a linear 
combination of i.y. Using conventional current density 
3u = iieh/2m)i^p,\/ij* - ^*V^P^) - ieym)A\i;,\^ with 
i = 1, 2, • • • iV, it is easy to show (e^//i) tu — J Ju ■ ds in 
the linear bias regime, where tu is the diagonal element 
of matrix T. In order to find an eigen-current density 3i 
such that 



-J- U^ I ii- ds 



(1) 



we define a generalized complex current density J^ so 
that (e^//i) tij = J ^ij ■ ds. The unitary transformation 
on the incoming wave function \Wa^m > suggests the 
following definition: 

where Tpi is the wave function in the scattering region 
due to the incoming state \Wa.i >■ From this defini- 
tion, we can prove the relationship (e^//i) tij = J 3ij ■ ds, 
as follows. It is not difficult to show that the following 
3i satisfies_Eq.([l|): J,_= {ieh/2m){4JiVi)* - ■4)*Vi}{) - 
(e2/m)A|V'jp where V'i = J2j^3Uji- Using this Ji, 
Eq.([T]) becomes 



2 /> 

jk 



(3) 



Denoting j7g the generalized current density matrix with 
matrix elements Jij , Eq. ([3]) is equivalent to 



e2 _ /• e^ 

— T^U^Ja-dsU or —T- 



Ja-ds. (4) 



We have also confirmed Eq. (|4]) numerically using specific 
examples including that shown in the inset of Fig[l}L. 
Therefore, to obtain eigen-current density matrix, we 
first diagonalize the transmission matrix T to find the 
unitary matrix U; we then calculate the generalized cur- 
rent density Jg according to Eq.®. The eigen-current 
density matrix is finally obtained by jTcigon = U^ JgU 
and plotted for visualization. 

Can eigen-channels be measured experimentally? To 
answer this question, as an example let's consider a two- 
probe device having two eigen-channels or two edge states 
in the presence of magnetic field. Assume one can per- 
form two experiments: (i) measurement of conductance 
and (ii) measurement of shot noise. Clearly, conduc- 
tance is given by G = ^(ti + ^2)- The shot noise is 

given by[i| S = j;;[ti{l - h) + i2(l - ^2)]- From these 
ti and ^2 can be determined. In the case of three eigen- 
channels, one needs to experimentally measure an addi- 
tional quantity [6], for instance the third cumulant of cur- 
rent Y =< AI{ti)AI{t2)AI{t3) >. In the linear regime, 
Y = (e^/h) X;, ti{l - U){1 - 2t,)X]- Hence by measuring 
G, S, and Y, one can determine ti,2,3- Therefore, the 
transmission eigen-channels are physical quantities mea- 
surable experimentally. 

Having prepared analytical tools, we now present nu- 
merical calculations on how the edge states are destroyed 
by increasing degrees of disorder. In the tight-binding 
representation, the Hamiltonian of 2D honeycomb lattice 
of graphene can be written as: 



H 



= Y. ^'^»"=« " * H 



^i2'K(f>^ 



'c\c 



i'^3 



(5) 



<y> 



where cj (ci) is the creation (annihilation) operator for an 
electron on site i. The first term in H is the on-site sin- 
gle particle energy where diagonal disorder is introduced 
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FIG. 1: (color online) (a,b) are "Phase diagram" of edge states 
of graphene for energies near the Fermi level (Dirac electrons): 
(a) in the {E, W) plane; (b) in the {B, W) plane, (c) The 
order of closing of eigen-channels at large disorder. Color 
coding are for number of eigen-channels. 



by drawing Ci randomly from a uniform distribution in 
the interval [—W/2, W/2] where W measures the disor- 
der strength. The second term in H is due to nearest 
neighbor hopping that includes the effect of a magnetic 
field. Here the phase 4>ij = J A-dl/(t)Q and 0o = h/e is 
the flux quanta. We fix gauge so that A = {By, 0, 0); and 
current flows in the x-direction. Transmission coefficient 
is given by T = Tr[T] where the transmission matrix T 
is obtained from T = TrCTlG" with G"'" being the 
retarded and advanced Green functions of the disordered 
scattering region. Quantities Fq {a = L, R) are the line 
width functions obtained by calculating self-energies S'' 
due to the semi-infinite leads [8|. The numerical data are 
mainly obtained from systems with 32x56 sites. In the 
calculations, energy and disorder strength are measured 
in unit of coupling strength t. 

Numerically, an edge state is identified if transmission 
coefficient of an eigen-channel is T > 0.999; if T < O.OOI, 
the eigen-channel is said to be closed. In addition, an 
edge state is said to be "destroyed" by disorder and be- 
comes a regular eigen-channel if its transmission T drops 
to below 0.999. Fig[T^ plots the phase diagram of edge 
states of graphene in the {E, W) plane with cj) = 0.0173 
and energy range [0, 0.85] where band dispersion is linear 
(Dirac electrons). A mesh of 600 x 480 points are scanned 
in {E, W) plane and up to 200 disorder configurations are 
averaged at each point. Several observations are in order. 
First, the edge states are destroyed one by one as W is in- 
creased. For instance, at E — 0.5 the v — 3 edge state is 
destroyed when W — 0.5. Very importantly, we empha- 
size that at this disorder, there are still three transmission 
eigen-channels although only two are edge states and the 
third being a regular eigen-channel having T < 0.999. In 
other words, the third eigen-channel is still there to par- 
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FIG. 2: (color online) Phase diagram of edge states in {E, W) 
plane for: (a) hole-like charge carriers in graphene; (b) elec- 
trons in a square lattice. The main difference between the 
two phase diagrams is due to opposite charges of hole and 
electron. 



FIG. 3: (color online) Current density of two edge states flow- 
ing from left to right along x-direction, for a square lattice of 
100 X 40 sites in the absence of disorder. Here E — 0.4 and 
<jf> = 0.052. Note that the color scale is different for the two 
panels. Insets: classical skipping orbits. 



ticipate transport although it is no longer an edge state. 
Increasing disorder to W^ = 0.7, the ly = 2 edge state is 
destroyed; finally when M^ = 1, all three edge states are 
destroyed. We note that edge states would be destroyed 
all at once if we had used the usual transmission coeffi- 
cient for each channel tu to characterize the edge states. 
Second, upon further increasing W, an insulating state 
is reached where all eigen- channels are closed. The order 
of channel closing is also one by one, in the same order as 
how edge states are destroyed. This is shown in Fig[TJ:. 
For instance, dX E — 0.5 there are three eigen-channels 
to start with, and at large disorder W — 3.1, one of them 
is closed leaving only two regular eigen-channels. Third, 
the edge states are easily destroyed at the subband edges 
while at the subband center they are most robust against 
disorder. This is because the energy of Landau levels is 
located at the subband edge. In the presence of disor- 
der, the Landau level is broadened with a finite width[l[. 
Hence the edge state that is close to one Landau level 
can easily relax toward it. Forth, it is more difficult to 
destroy an edge state at smaller energies. For Dirac elec- 
trons, the density of states is proportional to ^/E so that 
the level spacing of lower Landau levels is larger than the 
upper ones. For electrons with smaller energy it is far- 
ther away from nearby Landau level than electrons with 
larger energy. Hence a larger disorder is needed to relax 
the electrons to the nearby Landau level. Finally, Fig[T]3 
plots a "phase diagram" of edge states in the (B, W) 
plane for Dirac electrons. Once again, the edge states 
are destroyed one by one, similar to the phase diagram 
in the {E, W) plane. 

Fig[2^ depicts the "phase diagram" of edge states 
of graphene in (S, W) plane for higher energies in the 
range E = [2.545,2.97], where hole- like behavior occurs 
and band dispersion is non-linear (non-Dirac electrons). 
Again, edge states are destroyed one by one. While the 



"phase diagram" topology is similar to that for Dirac 
electrons (FiglT]), here the band dispersion is quadratic 
with equal energy spacing between the Landau levels. 
Due to this reason, the values of W that are needed 
to destroy the last edge state at different subband cen- 
ters are almost the same. We have also calculated the 
"phase diagram" of edge states of a square lattice in the 
(i?, W) plane using the same numerical method, results 
are shown in FiglJlD which are rather similar FiglJ^. In 
particular, it is more difficult to destroy edge states at 
low filling factor i/, consistent with the result of ReflSJ. 

Next, we examine the nature of those edge states that 
have survived disorder by calculating current density 
from Eqs.([2]) and ([4]) and plotting it along the propagat- 
ing direction (x-direction). Figl3] shows the current den- 
sity of two edge states in the absence of disorder for the 
square lattice model. Edge states are clearly seen. Since 
the two transmission eigen-channels have different longi- 
tudinal energies or effective velocities along the propagat- 
ing direction, it gives two different transmission patterns 
that correspond to two different skipping orbits of classi- 
cal trajectory[i3|. In FiglSK, current flows in the negative 
direction (blue region) near the sample boundary and in 
the positive direction (red region) away from it. There is 
a region between these opposite flows where the current 
density is very small. The classical trajectory of an elec- 
tron under Lorentz force is depicted in the inset, show- 
ing a nearly completed circular motion before colliding 
with the sample boundary. There is an one-to-one cor- 
respondence between the classical and quantum motion: 
near the sample boundary the flow is from right to left, 
while it flows opposite away from the boundary. Similar 
one-to-one correspondence is also seen in FiglSJD. For the 
same square lattice model, Fig|4]plots the current density 
of two eigen-channels for a particular disorder configura- 
tion W = \ where the eigen transmission coefficients are 
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FIG. 4: (color online) Current density along x-direction for 
the square lattice model as that of Fig|3] with disorder 
strength W = 1, energy E = 0.4 and magnetic field </> = 0.052. 
(a) For an edge state that has survived disorder; (b) edge state 
that has been destroyed by disorder. 



Ti = 0.9999 and T2 = 0.3385, respectively. In the numer- 
ical calculation, we have confirmed that the integral of J^ 
over any cross-section area along the propagating direc- 
tion gives the same value that is equal to {e^ /h)ti. From 
FigHK, it is obvious that Ti = 0.9999 is an edge state 
that survived this degree of disorder. Compare to FigU 
the pattern of current density with disorder scattering is 
clearly different. For the eigen-channel with T2 = 0.3385, 
it is clearly a non-edge state (FiglSb): there is a circu- 
lating patten with large current density in the middle of 
the scattering region, caused by the disorder scattering. 
Finally, we have also calculated current density for edge 
states in disordered graphene, and similar behaviors are 
observed as that of FigHl 

In summary, we have investigated the nature of edge 
states in disordered mesoscopic samples in the IQHE 
regime. Our results show that edge states are destroyed 
one by one as disorder strength is increased. In the insu- 
lating regime, all transmission eigen-channels are closed 
and the closing is also one by one in the same order as 
the edge states were destroyed. For graphene and the 
square lattice model, the "phase diagrams" have simi- 
lar topology but with some differences due to band dis- 
persions. We have introduced a quantity which is the 
generalized current density, using it the current density 



of each eigen-channel can be calculated in the presence 
of disorder, giving us a vivid physical picture on how 
edge states are destroyed. Since transmission coefficients 
of individual eigen-channels for mesoscopic samples can 
be determined experimentally in the IQHE regime — as 
we discussed in the paper, our conclusions on how edge 
states are destroyed by disorder should be testable ex- 
perimentally. 

We note in passing that how to generalize our results to 
the large sample limit is an interesting problem requiring 
further investigation. In particular, based on the pic- 
ture of Khmelnitskiill] and Laughlin[12], a global phase 
diagram of quantum Hall effect in the {B, W) plane was 
proposed by Kivelsonei.a?[13| for large samples which at- 
tracted considerable attention both theoretically!^, [3 
and experimentally |15j]. According to it, an integer quan- 
tum Hall state with a fixed filling factor v will floa t up 
in energy as the disorder strength increases [ill Il4l- I'^ 
the context of mesoscopic sample, this idea would mean 
the following. Consider Fermi level E — Ef, at the 
mesoscopic sample boundary there are, say, v edge states 
whose eigen- values cut this energy Ef. The "float up" 
idea means that when disorder is increased, the energies 
of edge states increase to higher values, i.e. they float 
up. Hence, at large enough disorder there will only be 
V — 1 edge states cutting Ef . This way, the system un- 
dergoes a series of transitions between i^ to i^ — 1 states 
etc.. Our results presented above, however, indicate that 
for mesoscopic samples edge states do not float up by dis- 
order, they are destroyed to become regular transmission 
eigen-channels which participate transport. 
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